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We explore a fluid description of the inflationary universe. In particular, we investigate a fluid 
model in which the equation of state (EoS) for a fluid includes bulk viscosity. We And that the three 
observables of inflationary cosmology: the spectral index of the curvature perturbations, the tensor- 
to-scalar ratio of the density perturbations, and the running of the spectral index, can be consistent 
with the recent Planck results. We also reconstruct the explicit EoS for a fluid from the spectral 
index of the curvature perturbations compatible with the Planck analysis. In the reconstructed 
models of a fluid, the tensor-to-scalar ratio of the density perturbations can satisfy the constraints 
obtained from the Planck satellite. The running of the spectral index can explain the Planck data. 
In addition, it is demonstrated that in the reconstructed models of a fluid, the graceful exit from 
inflation can be realized. Moreover, we show that the singular inflation can occur in a fluid model. 
Furthermore, we show that a fluid description of inflation can be equivalent to the description of 
inflation in terms of scalar field theories. 
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I. INTRODUCTION 

The various precise properties on inflation in the early universe [in has been revealed by the recent cosmological 
observations on the anisotropy of the cosmic microwave background (CMB) radiation such as the Planck satellite [j,@ 
and the BICEP2 experiment [ 1 , 0 ] j in addition to the Wilkinson Microwave anisotropy probe (WMAP) [ 1 , [3l|. The 
nature of inflation can be known from the spectrum of the primordial density perturbations |10l |. 

Recently, in Refs. [TTI.[T^. the convenient fluid description of the inflationary universe has been proposed in terms of 
scalar field theories, fluid models [ll|, and F{R) gravity theories [l^. Especially, the three observables of inflationary 
models, i.e., the spectral index of the curvature perturbations, the tensor-to-scalar ratio of the density perturbations, 
and the running of the spectral index, have been represented by using quantities in scalar field theorie^ fluid models, 
and EfR) gravity theories. Fluid models have been applied to cosmological issues such as inflation [isl - l^ and dark 
energy [ 24 ] (for recent reviews, see (l^. 

In this paper, by extending the preceding investigations of a fluid description [T^ for the inflationary universe, we 
derive the slow-roll conditions and construct the formulae of the observables for inflationary models of the spectral 
index Ug of the curvature perturbations, the tensor-to-scalar ratio r of the density perturbations, and the running as 
of the spectral index. We analyze the equation of state (EoS) for a fluid including bulk viscosity and examine a fluid 
model in which rig, r, and as can be compatible with the recent Planck results. 

In addition, we explicitly reconstruct the EoS for a fluid from the spectral index rig of the curvature perturbations. 
Particularly, we use the expression of Ug as a function of the number of e-folds N during inflation in the inflationary 
models including the Starobinsky inflation Q, from which the value of rig consistent the recent Planck analysis can 
be obtained. This reconstruction method for scalar field theories has been proposed in Ref. [2^. In this work, we 
present the reconstruction procedure in a fluid description. We also show that in the reconstructed models of a fluid, 
the slow-roll inflation, i.e., the de Sitter inflation, can occur. In these fluid models, the tensor-to-scalar ratio r of the 
density perturbations can meet the constraints acquired by the Planck satellite. The running ag of the spectral index 
can explain the Planck results. Moreover, we certify that in the reconstructed models of a fluid, the graceful exit from 
inflation can be realized. Furthermore, we demonstrate that in the fluid models, the singular inflation [25l - [^ can 
occur. In addition, we explain the equivalence between a fluid description of inflation and the description of inflation 
in terms of scalar held theories. We use units of fee = c = h = 1 and express the gravitational constant SttGn by 
with the Planck mass of Mpi = = 1.2 x 10^® GeV. 

The organization of the paper is the following. In Sec. II, we explain a fluid description of the inflationary universe. 
In particular, we present a fluid model in which the EoS for a fluid includes bulk viscosity and show that the fluid 
model can explain the recent Planck results of the three observables for inflationary models. In Sec. Ill, we reconstruct 
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the EoS for a fluid from the spectral index of the curvature perturbations. We certify that inflation can happen in the 
reconstructed models of a fluid, and that the tensor-to-scalar ratio of the density perturbations can be consistent with 
the Planck analysis. In Sec. IV, we investigate that the graceful exit from inflation can be realized in the fluid models 
reconstructed above. In Sec. V, we consider the singular inflation in a fluid model. In Sec. VI, we show that a fluid 
description of inflation can be equivalent to the description of inflation in terms of scalar field theories. Conclusions 
are presented in Sec. VII. In Appendix A, the slow-roll parameters in a fluid description are given. 


II. FLUID DESCRIPTION OF INFLATION 

We consider the case that the so-called slow-roll inflation driven by the potential V{(p) of a scalar filed (p occurs, 
which plays a roll of the inflaton field. We explain the procedure [llj | to represent the slow-roll parameters in terms 
of the Hubble parameter and its derivatives of the number of e-folds during inflation. Furthermore, with these 
representations of the slow-roll parameters, we describe the observables of inflationary models, namely, the spectral 
index of the curvature perturbations, the tensor-to-scalar ratio of the density perturbations, and the running of the 
spectral index in fluid models fl^ . 


A. Slow-roll parameters 


The action of 4> with the Einstein-Hilbert term is given by 

S = I -Vip)^ . 


(II.l) 


Here, g is the determinant of the metric and R is the scalar curvature. For the slow-roll inflation, the spectral 
index Ug of the curvature perturbations (i.e., the scalar mode of the density perturbations), the tensor-to-scalar ratio 
r of the density perturbations, and the running of the spectral index ag = dug/dlnk, where k is the absolute value of 
the wave number fe, are written as 


rig — 1 = — 6 e + 2g, r = 16e , ag 
where e, 77 , and ^ are the slow-roll parameters, defined as 

2k,^\ v{cj)) ) ’ K? ViP) ’ 


= 16e?7 - 24e2 - 2^^ ^ 


.2 _ 1 V'WV"'W 

{V{cP)f 


(H.2) 


(11.3) 


Here, the prime shows the derivative with respect to (p as V'{<p) = dV((p)/dp. Throughout this paper, the prime 
denotes the derivative with respect to the argument of the function, to which the prime operates. 

We take the flat Friedmann-Lemaitre-Robertson-Walker (FLRW) metric ds^ = —dt^ + a?{t) 2 3 {dx"^) i where 

a{t) is the scale factor. The Hubble parameter is defined hy H = d/a, where the dot means the time derivative. 

We express the slow-roll parameters in terms of H, which can be represented as H = H{N), namely, as a function 
of the number of e-folds N during inflation, defined as N = ln(af/ai) = J)/ Hdt, where Oi and Of are the values of 
the scale factor a at the initial time ti and the end time tf of inflation, respectively. To execute this task, with a new 
scalar field ip, we redefine p as p = p{ip), where (p is identified with N. We introduce a positive quantity ui{(p) (> 0) 
defined as uj{ip) = {dp/dip)'^, and represent V as a function of p, i.e., V{(p) = V {p{(p)). In the FLRW background, we 
derive the gravitational equations and rewrite them by using U!{(p) and V{ip). By solving the gravitational equations 
with respect to and V{tp), we obtain m 


Uj{'p) 


2 H'{N) 
h{N) 




Vip) 





W{N) 


N—ip 


(H.4) 


with H'{N) = dH{N)/dN. Here, the representations of = H{N) and (p = N are acquired as solutions for the 
equation of motion of p or p, and the gravitational field equations. It is seen from the first equation in (III.41) that 
since uj{p) > 0, we have H'{N) < 0. The slow-roll parameters in pi.31) can be rewritten by using uj{p) and Vipp). 
Accordingly, through the expressions of oj{p) and V{p) in (IH.4() . the slow-roll parameters can be described in terms 
of H{N) and its derivatives with respect to N. The resultant expressions have been given in Ref. EH- 
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B. Representation of a fluid 


For a general fluid model, the EoS is given by 


P{N) = -p{N)+f{p), 


(11.5) 


where p is the energy density of a fluid, P is its pressure, and f{p) is an arbitrary function of p. In the flat FLRW 
background, for such a fluid model, the gravitational field equations read 

^{H{N)f=p, ( 11 . 6 ) 

Ki 

-^HiN)H'iN) = p+P, (II.7) 

Since the EoS can be expressed as p{N) + P{N) = f{p), the second gravitational equation is rewritten to 
— (2/k^) H{N)H'{N) = f{p). Similarly, with the expression of the EoS shown above, the conservation law 
0 = p'{N) + 3{p{N) + P{N)) becomes 0 = p'{N) + 3/(p), where p'{N) = dp{N)/dN. From these second gravi¬ 
tational equation and conservation law, we acquire 
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V H{N) ) 


H"{N) 

H{N) 


3/'(p)/(p), 


(II. 8 ) 


with f'{p) = df{p)/dp. Owing to this equation, it is possible to express H{N) and its derivatives with respect to N 
only with p{N) and f{p{N)). Therefore, the slow-roll parameters can be described in terms of p{N) and f{p{N)), as 
is presented in Appendix A. As a result, by substituting the representations of the slow-roll parameters in Appendix 
A into the expressions of observables of the inflationary models in (III.31) . we obtain the fluid description of rig, r, and 
Ofs- In Ref. [^, the explicit expressions of rig, r, and ofg have been shown^. 

The form of the EoS for a fluid can also be represented as w{N) = P{N)/p{N) = —1 -|- f{p)/p{N), from which we 
find f{p)/p{N) = w{N) + 1. When \f(p)/p(N)\ 1, and f{p) and p varies very slowly in the inflationary stage, the 

approximate expressions of Ug, r, and Og read [l2l| 


s) (1-6 


fip) r,, fip) n f f(p) 


24 


-9 


p{N) ’ p(N) ’ yp{N) 


= (1 - 6 {w{N) + 1) , 24 (w(N) + 1) ,-9 {w(N) + 1)^), 


(II.9) 

(II.IO) 


where in deriving (III.101) . we have used the relation f{p)/p{N) = 'w{N) -|- 1. 


C. Fluid model in which the EoS for a fluid includes bulk viscosity 

We investigate a fluid with the following EoS 

P = -P + Ap^ + C{H), (II.ll) 

where A and j3 are constants, and C{H) is bulk viscosity. As a specific case, we consider that (^{H) has the following 
form 


an) = , ( 11 . 12 ) 

where ( and 7 are constants. We note that the mass dimension of A is —4 {fd — 1), whereas that of C is — (7 — 4). 
From the Friedmann equation (III.61) for the expanding universe {H >0), we get H = (k/ Vs) ^/p. Hence, C,{H) can 
be written as a function of p, namely, C{H) = ({H{p)). Consequently, by comparing Eq. (III.5() with Eq. (Ill.llj) and 
using Eq. (III.12|) . we acquire 

f{p) = Apf^ + aH{p)) =Ap^ + C ^ . (11.13) 


^ The other way to describe oig has been examined in Ref. [^ . 
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Here, we state a physical reason why we have considered the case that C(^) is expressed by a power in H as 
given in Eq. (III.121) and hence f{p) is represented by the linear combination of two kinds of a power in p as shown 
in Eq. (III. 131) . It is considered that only for such a case that f{p) is expressed by a series of a power in p, through 
a phenomenological approach, it is possible to analytically study the quantitative features of the EoS for a fluid to 
realize inflation in which the three observables of inflationary models, namely, the spectral index of the curvature 
perturbations, the tensor-to-scalar ratio of the density perturbations, and the running of the spectral index, can 
explain the recent Planck results, as is demonstrated below. 

The Planck analysis [1,0 has shown that rig = 0.968 ± 0.006 ( 68 % CL), r < 0.11 (95% CL), and ag = —0.003 ± 
0.007 ( 68 % CL). If f{p)/p{N) = 4.35 x 10“^, i.e., w = —0.996, from Eq. (111.91) or Eq. (III.101) . we have (ng,r, ag) = 
(0.974,0.104,-1.70 x 10“"^). These results are consistent with the Planck data. 

We here demonstrate that it is possible to realize these Planck results by choosing appropriate values of the model 
parameters in Eq. (III.Ill) with Eq. (111.121) . In other words, we explicitly derive the values of the model parameters 
leading to f{p)/p = 4.35 x 10“^. It follows from Eq. (111.131) that 


f{p) 

P 


M 


fi-i / P_ 


+ c 


(V 3 


n^/2-l 


7 / 2-1 


Ap: 




ginf \ 

Ho J 




\V3j 


Ho ) 


(11.14) 

(11.15) 


where Pc = = 8.10 x 10“'*^GeV^ is the critical density. Ho = 100hkmsec“^ Mpc“^ = 2.13/i x 10“'^^GeV with 

h = 0.678 [1] is the current Hubble parameter [s^, and iLinf is the Hubble parameter at the inflationary stage. For 
simplicity, we set 7 = 2/3. In this case, from Eq. (III.151) . we obtain 


fip) 


= J 


H 


inf 


J = 


Ho 

A + C 


2(/3-l) 


(^) 


2/3' 


3-1 


(11.16) 

(11.17) 


For the simplest case that /3 = 1, when J = 4.35 x 10“^, regardless of the scale of inflation iLinf, the Planck results 
can be realized. Moreover, in the case that /3 = 2, for example, if (iLinf, J) = (1.0 x 10^°GeV, 9.10 x 10“^°^), 
(1.0 X 10® GeV, 9.10 x 10“®^), we can explain the Planck data. 


III. RECONSTRUCTION OF THE EOS FOR A FLUID FROM THE SPECTRAL INDEX 

In this section, we reconstruct the EoS for a fluid from the spectral index of the curvature perturbations. Such a 
reconstruction has been studied for the case of scalar field theories in Ref. [^ . 


A. Reconstruction procedure in a fluid description 


For the slow-roll inflation in scalar field theories, whose action is given by Eq. (IH.1|) . the spectral index rig of the 
curvature perturbations, the tensor-to-scalar ratio r of the density perturbations, and the running ag of the spectral 
index are is derived as follows [ 2 ^: 


rig 



In 


( 1 dViN) \- 

\V^{N) dN )\ ’ 


8 dV{N) f ^ dV{N)\ 


(HLl) 


Similarly to the case of scalar field theories, in a fluid model, it is possible to reconstruct the EoS for a fluid from 
the spectral index rig of the curvature perturbations. If we have the form of rig as a function of N, by using the 
first relation in (IIII.ll) , we can obtain the expression of V(N). Thanks to the Friedmann equation (III.61) , the Hubble 
parameter is related to V{N), and hence we get H = H{N). In a fluid model, with the other gravitational field 
equation (III.71) . we can acquire the form of f{p) through the EoS in Eq. (III.51) . 
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B. Inflationary models with rig — 1 = —2/N 


We demonstrate the reconstruction procedure in a fluid description by exploring the inflationary models in which 
rig is given by 


n, — 1 = — — . 

N 


(IIL2) 


It is known that in the Starobinsky iirflation {B? inflation) Q, rig and r are expressed as the relation (IIII.2|) and 
r = respectively. If = 60, we And Ug = 0.967 and r = 3.33 x 10“^, which are corrsistent with the Planck 

data 0 (for a recent detailed review of inflation in modified gravity theories, see, for instance, [l^). The relation 
pil.2l) can be satisfied also in the chaotic inflation and the Higgs inflation with its non-minimal gravitational 
coupling [s^, or the so-called a-attractor [s^, which connects the Starobinsky, quadratic chaotic, and Higgs inflations. 
By combining the relation (IHI.2I) with the first equation in (IHI.II) . we And 


V{N) = 


(HI.3) 


{Ci/N) + C2 ’ 

with Ci(> 0) and C 2 constants, the mass dimension of which is four. For the potential V{N) in Eq. (IHI.3I) . from the 
second relation in piI.lL the tensor-to-scalar ratio r of the density perturbations is expressed as 

(HI.4) 


lV[l + (C2/Ci)iV] ’ 

Furthermore, with the third relation in (IHI.IIl . the running Og of the spectral index is written as 

2 


Qg - 


1V2 


(HI.5) 


By using this expression, for N = 60, we acquire Og = —5.56 x 10 This value is consistent with the Planck analysis. 

In a fluid model, instead of the inflaton potential F, we use the EoS in Eq. (III.51) . In the FLRW background, the 
Friedmann equation (III.61) is written as {^/k^) iH{N))'^ = p{N) m V{N), where the last approximate equation follows 
from the slow-roll approximation that the kinetic term is much smaller than the potential one as ( 1 / 2 ) ^2 y 
From the relation p K.V with Eq. (IHI.3I) , we have 

Furthermore, it follows from the Friedmann equation with the slow-roll approximation shown above that the Hubble 
parameter is expressed as 


H{N) 


1 


y 3[(Ci/iv) + c2] ’ 

where {Ci/N) + C 2 > 0. From Eq. (|H. 6 I) and (III.71) with the Hubble parameter in Eq. (IHI.7I) . we obtain 

P = -p- ^HiN)H\N) ^-p- 


By comparing Eq. (III.51) with Eq. (IHI. 8 |) . we acquire 

^ ~ 4 

M - 




-^(1-2C2P + C2V) . 


(HI.7) 


(HI. 8 ) 


(HI.9) 


Here, in deriving the second approximate equality, we have used the Friedmann equation pi. 6 |) and Eq. (IHI. 6 p . 


C. Fluid models and inflation 

Next, we explicitly show the models of a fluid, in which the values of rig and r are consistent with the Planck results. 
Plugging Eqs. dlLTTI) and (limi) . we have the form of EoS for a fluid 

P = -P + f{p) = -p + V + C • (ni.io) 

With the results in the preceding subsection, we decide the models parameters: A, /3, and 7 , in which the relation 

(IHI.2I) can be satisfied. 
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1. Case (i): \C 2 p\ S> 1 


When \C 2 p\ ^ 1, from Eq. (IIII.9I) . we have 


Up) 


2C2 2 

3Ci ^ sc/ ■ 


(III.ll) 


Since the value of N given by Eq. (Illl.hj) has to be positive, we find C 2 < 0- In addition, 
during inflation has to be much larger than unity such as iV = 60, and hence, from Eq. 
\C 2 p\ ^ 1, we acquire (—^ 2 ) /Ci k. 1/N <C 1. From Eqs. (IIII.IOI) and (IIII.llI) . we get 


the number of e-folds N 
pil. 6 |) and the condition 


where in deriving the second approximate equality, we have used (— C'2) /Ci « 1/iV. For example, if |C'2p| = 0(10) 
and (—C'2) /Cl ~ 1/N, where, e.g., N > 60, from Eq. pII.12p . we can obtain w ^ —1. This implies that the slow-roll 
inflation, namely, the de Sitter inflation, can occur, and hence the scale factor can be represented as 


a{t) = Oi exp [i^inf (t — ti)] . (III.13) 

Moreover, if (—C 2 ) /Ci < 1/fV, from Eq. (IIII.4|) . it is seen that for N > 73, the tensor-to-scalar ratio r of the density 
perturbations can meet r < 0.11, which is consistent with the Planck results. 

Through the comparison between this expression and Eq. (III.131) . we see that these expressions become equivalent, 
i.e., the linear combination of p and p^. In this case, there are two combinations of the model parameters, which will 
be called as Model (a) and Model (b) as follows 

‘ic' - 

Model (a): A=^, C =- 7 ^^, /3 = 1, 7 = 4, (III.14) 

oOl Oi K, 

and 

(^2 _ 2(^ 

Model (b): ^ = C = 7^4 ’ /3 = 2, 7 = 2. (III.15) 

oOl C'l/'C 

In piI.14(l and piI.15(l . when the second relations have been derived by using the forth relations. As a result, the 
EoS of a fluid can explicitly be reconstructed. 


2. Case (ii): |C 2 p| <C 1 


On the other hand, if |C 2 p| -C I, it follows from Eq. (IIII.Qj) that 


f{p) 


I 2 C 2 

3Ci 3Ci ^ 


(III.16) 


With Eq. (IIII. 6 I) and the condition |C 2 p| <C I, we have Cipfv N I, and eventually we also find IC 2 I /Ci <C I. From 
Eqs. (IIII.IOI) and (IIII.llI) . we acquire 


P 


w = — 


P 


1 1 

3^ 


2 

3 




(III.17) 


Here, the second approximate equality follows from Cip ~ N. Accordingly, from Eq. (IIII.17I) with 1/N <C 1 and 
\C 2 \/C\ -C 1, we see that w si — 1 can be met. Thus, the slow-roll (de Sitter) inflation can happen, and the scale 
factor can be expressed by Eq. (IIII.13|) . In addition, for C '2 > 0 and C 2 /C 1 < 1/N, by using Eq. (IIII.4I) . it is found 
that even for N > 60, the tensor-to-scalar ratio r of the density perturbations becomes r < O.ll, which is consistent 
with the Planck results. On the other hand, for C 2 < 0 and IC 2 I/C 1 < 1/A^, it follows from Eq. (IIII.4I) that for 
N > 73, we can get r < 0.11, similarly to that in Case (i) described above. 

The comparison of this expression with Eq. (III.131) leads to the following combinations of the model parameters, 
which will be named Model (c) and Model (d) as follows 


2 C 2 


Model (c) : A 


1 


Cik2 ’ 


/3 = 0 , 7 = 2 


(III.18) 
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Table I: Fluid models with the EoS in Eq. (IIII.IOI) realizing the relation (IIII.2I) . Here, C\ > 0. In Case (i), \C 2 p\ 3> 1 and 
C 2 < 0, whereas for Case (ii), \C 2 p\ <S 1 and C 2 can take both the positive and negative values. 


Case 

Model 

A 

c 

P 

7 

(i) 

(a) 

2(72/ (3(7i) 

-3(71/ {Cik'^) 

1 

4 

(i) 

(b) 

-G|/(3Ci) 

2(72/ 

2 

2 

(ii) 

(c) 

-1/ (3Ci) 

2(72/ 

0 

2 

(ii) 

(d) 

2(72/ (3(7i) 

-l/(3Ci) 

1 

0 


and 


IC - 1 

Model (d): ^ = 

In Eqs. (IIII.18I) and (IIII.19I) . with the forth relations, the second ones have been derived. In Table HI the fluid models 
with the EoS in Eq. (IIII.IOI) satisfying the relation (IIII.2I) are summarized. 

We remark that if C 2 >0, the inflaton potential can correspond to the one in the Starobinsky inflation. From the 
investigations in the scalar field theories, we have C 2 = (2/3) C\ [IJ. In this case, for the models in Eqs. (IIII.18I) and 
(IIII.19I) . we obtain (/ = 4/ (3 k^) and A = 4/9, respectively. 


IV. GRACEFUL EXIT FROM INFLATION 

In this section, we examine whether the graceful exit from inflation can occur in a fluid model. We analyze the 
instability of the de Sitter solution {H = Hmf (> 0) = constant) during inflation by taking the perturbations of the 
Hubble parameter as follows 


H = i7i„f + iJinf(5(t) . 


(IV.l) 


Here, |(5(t)| 1, and hence HiniS{t) denotes the perturbations from the de Sitter solution Hinf- 

We rewrite Eq. (jH.8ll as the following second differential equation with respect to the cosmic time t: 











= 0 . 


We define the form of S{t) as 


(IV.2) 


Sit) = e^*, 


(IV.3) 


where A is a constant, so that we can investigate the instability of the de Sitter solution. If there is a positive solution 
of A, the de Sitter solution can be unstable. Therefore, the universe can exit from inflation, and the reheating stage 
can follow, because the absolute value of d(t) with A > 0 becomes larger as the cosmic time grows at the inflationary 
stage. 

We substitute Eq. (IIV.II) with Eq. (IIV.3I) into Eq. (IIV.2I) and take the first order of 6(1). Accordingly, we get 


A^ 


1 




2 = 0 , 


2 = ^ (4/3 - 1) A^ 



(/3+I) (2/3 + 7-1)AC 


(IV.4) 

+ i (27 - 1) ■ (IV.5) 


We see that the solutions of Eq. (IIV.4I) are given by 


A = A± = i 


1 

Hinf 


72 . 


(IV.6) 


If 2 > 0, we can acquire the positive solution of A = A+ > 0. As a result, the exit from inflation can gracefully occur. 





































Table II: The EoS for the fluid models reconstructed in Sec. Ill and the conditions that in these models, the graceful exit from 
inflation can be realized. In these models, ris — 1 = —2/N = 0.967 for N = 60, r < 0.11 for > 73 in Models (a), (b) and 
Models (c) and (d) with C 2 < 0, or N > 60 in Models (c) and (d) with C 2 > 0, and as = —2fN^ = —5.56 x 10“^ for = 60 
can be realized. These values can explain the Planck data. Legend is the same as Table H] 


Case Model 

EoS 

Conditions for the graceful exit from inflation 

(i) (a) 

P=-p+ [ 2 C 2 / (3Ci)] p - \‘iCll (CiK^)] 

No condition 

(i) (b) 

P = -p- [Ci/ (3Ci)] p2 -k [2(72/ (Ci«2)] 

No condition 

(ii) (c) 

P=-p- [1/ (3Ci)] -t [ 2 C 2 / (CiK^)] 

C 2 < 0 or C 2 > (1/36) (K/Hintf 

(ii) (d) 

P=-p+[2C2/(3Ci)]p-[l/(3Ci)] 

C 2 < 0 or C 2 > (1/18) («/Pinf)" 


Concretely, in the fluid models reconstructed above and summarized in Table HI we check whether the graceful exit 
from inflation can be realized or not, namely, whether Q can take a positive value or not. If the universe cannot 
successfully exit from inflation, inflation does not ends, and therefore such a scenario corresponds to the so-called 
eternal inflation. By substituting the values of A, /3, and 7 in Models (a), (b), (c), and (d), given by Eqs. pil.141) . 
(IIII.15|) . (IIII.18I) . and (IIII.19I) . respectively, into Eq. (IIV. 6 |) . we obtain the expressions of Q in each models. To evaluate 
the values of Q, we take into account the following facts. For all of the models, Ci > 0. On the other hand, in Models 
(a) and (b), C 2 < 0, while in Models (c) and (d), C 2 can become both the positive and negative values. For Models 
(a) and (b) in Case (i), we And 


Model (a) : Q = 2 


C 2 


CiJ V 


/ gjnf 


6 - 45C2 




Model (b) : Q = 6 ( ^ 


Hinf \ 


2 - 15^2 




21^1 , 


/ gjnf 


> 0 , 


> 0 . 


(IV.7) 

(IV. 8 ) 


Accordingly, we always have Q > 0. While, for Models (c) and (d) in Case (ii), we acquire 


Model (c) : Q = 


Co 


CiJ V 


/ gjnf 


Model (d) : Q 


=2(^y(^ 


Cj \ K 


3 C 2 

6 


12 


i7i„f 


ginfV 1 
K J 3C2 


(IV.9) 

(IV.IO) 


From these relations, we find that if C 2 < 0, we get Q > 0, whereas, in the case that C 2 > 0, if the following conditions 
are satisfied 


C 2 > 


C 2 > 


1 f K 


36 Vi7inf 
1 f K 


18 VFIinf 


for Model (c), 

(IV.ll) 

for Model (d), 

(IV.12) 


we can obtain Q > 0. Thus, for the reconstructed models of a fluid in the previous section, it is possible for the 
universe to gracefully exit from inflation. 

In table El we present the summary of the reconstructed fluid models. We show the EoS of these models in the 
form of Eq. (III.Ill) so that a term inspired by bulk viscosity can clearly be seen. In these models, the three observables 
of inflationary cosmology can be compatible with the Planck results. First, the spectral index Us of the curvature 
perturbations is expressed as rig — 1 = —2/N in Eq. (IIII.2I) . which can lead to 0.967 for N = 60. Second, The 
tensor-to-scalar ratio r of the density perturbations can satisfy the upper limit of r < 0.11. In Models (a) and (b) 
[Case (i)] and Models (c) and (d) [Case (ii)] with (72 < 0, if Al > 73, we can obtain r < 0.11. While, in Models (c) 
and (d) [Case (ii)] with C 2 > 0, when N > 60, we can find r < 0.11. Third, the running ag of the spectral index is 
given by ag = —2/N'^ in Eq. (IIII.5|) . From this expression, we have ag = —5.56 x 10“^. These values of rig, r, and 
Qfg are consistent with the Planck results. Moreover, the universe can gracefully exit from inflation. We describe the 
conditions for the graceful exit from inflation. 
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V. SINGULAR INFLATION IN A FLUID MODEL 

In this section, we study the singular inflation in a fluid model. In this inflationary scenario, the idea of 
finite-time future singularities in the context of the dark energy problem is applied to inflation in the early universe. 

The finite-time future singularities are classihed into four types [s^ . Their features in modihed gravity theories 
have also been analyzed in detail (for a detailed review on the finite-time future singularities, see [l^)- Among 
them, the formulation of the Type IV singularity can be used in the singular inflation, because there is no divergence 
in terms of the scale factor, the energy density and pressure of the universe. 

In the Type IV singularity, for t —>■ tg, where tg is the time when the singularity appears, the scale factor a, the 
effective (i.e., total) energy density pes and pressure PeS of the universe become finite as a —>■ Og, pes —t 0 and 
iT’effI —>■ 0. Here, Og is the value of a at t = tg. The case that peS and/or |Peff| become non-zero finite values at 
t = ts [1^ is also included in the Type IV singularity. However, the higher derivatives of H diverge. 

In the following, we explore the inflationary stage in which there is only a component of a fluid. Therefore, for 
simplicity, we describe peS and Peff by p and P, respectively. We consider the case that the Hubble parameter and 
scale factor during inflation are represented as 


H = Pinf + H (4 


a = a exp 


Hi-ait 


H 


q>l, 

{u - 


(V.l) 

(V.2) 


where P, q, and a are constants, and the mass dimension of P is q+ 

In the flat FLRW universe, from the gravitational field equations, the energy density and pressure of the universe 
are given by 


P = 


3P2 


P = - 


2P -h 3P2 
^2 


(V.3) 


It is seen from Eq. (jV.2ll and the expressions in (IV.31) with Eq. (IV.Ill that in the limit t ^ ts, all of a, p, and 
P asymptotically approach finite values, while the higher derivatives of P diverge. Thus, the Type IV singularity 
appears at t = ts- By using the expressions of p and P in (IV.21) with Eq. (IV.II) . we find the following EoS for a fluid 


P = -P + f{p). f{p) 




(<?-!)/<? 


(V.4) 


Here, f{p) can be described as the series of power of p. In fact, if Pinf/ \/tP’pji = Pinf/P ^ 1, where the equality 
comes from the first equation in (IV.31) . we find 


/(P) 


2 pi/« 

3('?-i)/(2'j) t^<,q+^)/q 


p(9-i)/(29) 


■\/3 (g 1) Pjnf ^-i/( 2 g) 
q K 


(V.5) 


where we have taken the first order of the quantity 


. It follows from Eq. (IV.51) that f{p) is represented 


as a linear combination of two powers of p, similarly to that in Eq. (jH.13() or Eq. (jlll.lip . Hence, this model can be 
regarded as a kind of the fluid models reconstructed in Sec. HI. 

Furthermore, Eq. (IV.51) divided by p reads 


fjp) 

P 




^ (g - 1) Pinf 

V 3 ) 

q ^Jk^pIZ_ 


3 


f H y/" 

r (g-l)Pinf] 


[ q H \ 


(V.6) 


Here, in deriving the last equality, we have used the first equation in (IV.31) . It is seen from Eq. (IV.6p that for 
P/P'J+^ 1, we have f{p)/p ^ 1. In this case, the observables of the inflationary models, i.e., rig, r, and Og, 

can approximately be represented by Eq. (III.91) . and the values of Ug, r, and Og can be compatible with the Planck 
analysis, as stated in Sec. H C. 

We explain the existence of limit of C = 0 in Eq. (III.121) . in which the term C{H) in Eq. (IH.12I) will not exist and 
therefore the EoS for a fluid in Eq. (III.Ill) reads P = —p + Ap^. In such a limit, from Eq. (III.51) . we have f{p) = Ap^, 
namely, the term f{p) consists of the single power of p. On the other hand, in Eqs. (|V.5p and (jV.bp . the form of f{p) 
is a linear combination of two kinds of power of p. The form of f{p) can approximately be given in Eqs. (IV.51) and 
(IV.61) only if the singular inflation occurs and the Hubble parameter and the scale factor are expressed as Eqs. (IV.II) 
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and (IV.21) . respectively, This means that for a fluid without the term C{H) in Eq. (III.121) . the singular inflation cannot 
be realized. Thus, the existence of the term CiH) can influence the dynamics of the universe filled with a fluid in the 
early universe. 

As a consequence, it is considered that the singular inflation can be realized in the fluid models in which the spectral 
index of the curvature perturbations can explain the recent Planck results. 


VI. EQUIVALENCE BETWEEN A FLUID DESCRIPTION OF INFLATION AND THE DESCRIPTION 
OF INFLATION IN TERMS OF SCALAR FIELD THEORIES 


In this section, we demonstrate that a fluid description of inflation can be equivalent to the description of inflation 
in terms of scalar field theories (for further related investigations, see Ref. [13 )■ The action of scalar field theories is 
expressed as 


5 " = 





(VI.l) 


Here, u}{^) is a coefficient function of kinetic term of the scalar field tp and V{(p) is the potential of p>. Starting from 
a fluid description, we construct a scalar field theory with the same EoS as that of a fluid. By this process, we obtain 
the expressions of and V{ip) of the corresponding scalar field theory to a fluid description. Consequently, we can 
represent a fluid description as the description of a scalar field theory. 

It is known that in the FLRW background, uj{ip) and Vij-p) can be described as [40j | 


Ul{ip) 

V{p}) 


2 dJ{ip) 
dp ’ 


dJ{p) 

dp 


(VI.2) 

(VI.3) 


with J(pp) an arbitrary function of p. Here, we can take p = t and H = J(t) because p can be treated as an auxiliary 
scalar quantity. On the other hand, the energy density p and pressure P of the scalar field p read 


P = ^w((p)(p2 + V((p), 

(VI.4) 

p = 

(VI.5) 

With these equations, we find that oj{p) and V{p) are given by 

ujp^ = p + P = f{p ), 

(VI.6) 


(VI.7) 


In deriving the second equalities in Eqs. (IVI.6I) and (I VI. 71) . we have used the EoS of P = —p+f{p) in a fluid description 
in Eq. (III.51) . From the Friedmann equation (III.61) . we find p = 2>H^/. Therefore, when we have H{= we can 

express p = p{t{p)) = p{p) as a function of t{= p). Eventually, from Eqs. (IVI.6I) and (IVI.7I) with p = p(p), we can 
acquire the expressions of w = u}{p) and V = V(p). By using these processes, we can obtain the description of a 
scalar field theory corresponding to an original fluid description. 

Moreover, we consider the opposite approach from the description of a scalar field to a fluid description. We first 
have a scalar field action with uj{p) and V{p) in Eq. (IVI.ll) . It follows from Eqs. (I VI. 41) and (I VI. 51) with 4> = t and 
H = J{t) that the EoS of the universe w = P/p = —1 + f{p)/p in Eq. (jVI.ll) with Eq. (III.51) . Plugging this relation 
with Eq. (IVI.4I) and comparing the relation obtained with re = —1 + /(p)/p, we get f{p) in a fluid description. Thus, 
it can be considered that both approaches shown above suggests the equivalence between a fluid description and the 
description in terms of scalar field theories. 

As a fluid description, we the following case that f{p) is given by 


f{p) = fi + h , (VI.8) 

where /i, / 2 , and u are constants, and p* is a fiducial value of p For /(p) in Eq. (IVI.8I) with u = 1, we have the 
following EoS 




(VI.9) 
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By comparing this form with Model (d) in Table HIl it is seen that if / 2 /p* = 2(72/ (3(7i) and fi = — [1/ (3(7i)], the 
form of the EoS in Eq. (I VI. 91) is equal to Model (d). This means that in a fluid description with the form of the EoS in 
Eq. (IVL8(I . the spectral index of the curvature perturbations, the tensor-to-scalar ratio of the density perturbations, 
and the running of the spectral index can be consistent with the Planck results. 

From the substitution of Eq. (IVI.8|) into the relation re = —1 + /(p)/p, Eqs. (IVI.6I) and (IVI.7|) . we acquire 


w 



P 


!i 





fl + h 



V = p 


1 

2 


/ 1+/2 



(VI.IO) 

(VI.ll) 

(VI.12) 


In addition, for instance, we consider the case that the Hubble parameter and the scale factor during inflation are 
expressed by 


H (VI.13) 

a = at^, (VI.14) 


where 1 and d(^ 0) are constants. Such a case of » 1 corresponds to the quasi-de Sitter inflation (i.e., the 
slow-roll inflation). In this case, it follows from Eqs. (I VI. 21) and (IVI.3I) with (p = t and H = J{t) that 


w(<p) 

V{p) 


2h 1 

h (3h — 1 ) 1 
^2 


(VI.15) 

(VI.16) 


As described above, when the Hubble parameter H can be represented as a function of t, from Eqs. (IVI.2|) and (IVI.3I) 
with if = t and H = J(t), the expressions of a; = w((p) and V = H(p) can be derived explicitly. 


VII. CONCLUSIONS 

In the present paper, we have investigated the description of the inflationary universe in the framework of a fluid 
model in which the EoS for a fluid includes bulk viscosity. It has been found that in a fluid description, the three 
observables of inflationary models, namely, the spectral index rig of the curvature perturbations, the tensor-to-scalar 
ratio r of the density perturbations, and the running Og of the spectral index, can be consistent with the recent Planck 
results. 

Furthermore, we have explicitly reconstructed the EoS of a fluid model from the spectral index rig of the curvature 
perturbations. Particularly, we have used the expression of rig as a function of the number of e-folds JV in the 
inflationary models, where the value of rig can explain the Planck data, including the Starobinsky inflation. It has 
been shown that for the fluid models reconstructed from the spectral index, indeed, the slow-roll (de Sitter) inflation 
can occur. It has also been certified that in these fluid models, the tensor-to-scalar ratio r of the density perturbations 
can meet the upper limit found by the Planck analysis. The running Og of the spectral index can be compatible with 
the Planck results. 

In our previous work [I^, since we have considered a fluid without the term ((ff) in Eq. (III.121) . only for the special 
case that the EoS for a fluid is approximately equal to — 1 as w = P/p k, —I, it has been shown that in a fluid model, 
the three observables of inflationary models can be consistent with the Planck results. On the other hand, in this 
work, we have introduced the additional term C,{H) in Eq. (III.121) into the EoS for a fluid as in Eq. (III.111) . As a 
result, it has been found that also for cases in which the value of the EoS w for a fluid is apart from —1, in such a 
fluid model, the three observables of inflationary models can be compatible with the Planck analysis. 

In addition, we have examined the instability of the de Sitter solution at the inflationary stage by analyzing the 
perturbations of the Hubble parameter. It has been performed that the universe can gracefully exit from inflation 
in the reconstructed models of a fluid. We have also derived the conditions for the graceful exit from inflation to be 
realized in the reconstructed fluid models. 

Moreover, we have explored the singular inflation in a fluid model by using the formulations to describe the Type 
IV singularity, which is one of the four types of the finite-time future singularity. It has been demonstrated that the 
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singular inflation can be realized in the fluid models where the spectral index of the curvature perturbations can be 
compatible with the Planck data. 

It has also been studied that a fluid description of inflation can be equivalent to the description of inflation in terms 
of scalar field theories. 

Consequently, not only the representation of inflation in scalar field theories but also a fluid description of the 
inflationary universe can explain the observational results acquired by the Planck satellite. 

The present method of the reconstruction may equally be applied in the case that the universe is hlled with several 
coupled fluids. This description may also be applied to the cosmological evolution from modified gravity consistent 
with a fluid description at the background evolution level. 
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Appendix A: Slow-roll parameters in a fluid description 

In this Appendix, we present the slow-roll parameters in a fluid description. The explicit expressions are given by 

e = 7:P[N)f{p) 


2p{N) - f{p) 


_ 3p{N) I , _ 5 f{p)f{p) ff{p)Y 

2p(iV)-/(p) + 2 p{N) 


1 p'jN) 

3 f{p) 


{ I \ ^ r>/*(p)/ (P) I ( f i.P') 

(/(p)) +f(p}f + 
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-f 
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2i2piN)-fip))^ 
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Ykw) V - 
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